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Spin Models on Thin Graphs 

C. F. Baillie a and D. A. Johnston 13 

a Computer Science Department, 
Campus Box 430, University of Colorado, 
Boulder CO80309, USA 

b Mathematics Department, Heriot-Watt University, 
Edinburgh, EH14 4AS, United Kingdom 

We discuss the utility of analytical and numerical investigation of spin models, in particular spin glasses, on 
ordinary "thin" random graphs (in effect Feynman diagrams) using methods borrowed from the "fat" graphs of 
two dimensional gravity. We highlight the similarity with Bethe lattice calculations and the advantages of the 
thin graph approach both analytically and numerically for investigating mean field results. 



1. INTRODUCTION 

The analytical investigation of spin glasses on 
random graphs of various sorts has a long and 
honourable history though there has been 

little in the way of numerical simulations. Ran- 
dom graphs with a fixed or fixed average con- 
nectivity have a locally tree like structure, which 
means that loops in the graph are predominantly 
large, so Bethe-lattice-like || (ie mean field) criti- 
cal behaviour is expected for spin models on such 
lattices. Given this, the analytical solution for 
a spin model or, in particular, a spin glass on 
a Bethe lattice |||| can be translated across to 
the appropriate fixed connectivity random lattice. 
Alternatively, a replica calculation can be carried 
out directly in some cases for spin glasses on var- 
ious sorts of random lattices. 

A rather different way of looking at the problem 
of spin models on random graphs was put forward 
in ||, where it was observed that the requisite 
ensemble of random graphs could be generated 
by considering the Feynman diagram expansion 
for the partition function of the model. For an 
Ising ferromagnet with Hamiltonian 



<ij> 



(1) 



where the sum is over nearest neighbours on 
three-regular random graphs (ie (f> 3 Feynman di- 



agrams), the partition function is given by 

where N n is the number of undecorated graphs 
with 2n points, K is defined by 
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and the action itself is 



a,b 
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(2) 



(3) 



where the sum runs over ± indices. The coupling 
in the above is g = exp(2/3J) where J = 1 for 
the ferromagnet and the <p+ field can be thought 
of as representing "up" spins with the </>_ field 
representing "down" spins. An ensemble of z- 
regular random graphs would simply require re- 
placing the </> 3 terms with <fi z and a fixed average 
connectivity could also be implemented with the 
appropriate choice of potential. 

This approach was inspired by the considerable 
amount of work that has been done in recent years 
on N x N matrix [] versions of such integrals which 
generate "fat" or ribbon graphs graphs with suf- 
ficient structure to carry out a topological expan- 
sion |Q because of the matrix index structure. 

N, the size of the matrix is not to be confused with n, 
the number of vertices in the graph! 
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The natural interpretation of such fat graphs as 
the duals of triangulations, quadrangulations etc. 
of surfaces has led to much interesting work in 
string theory and particle physics || . The parti- 
tion function here is a poor, "thin" (no indices, 
so no ribbons), scalar cousin of these, lacking the 
structure to give a surface interpretation to the 
graph. Such scalar integrals have been used in 
the past to extract the large n behaviour of vari- 
ous field theories Q again essentially as a means 
of generating the appropriate Feynman diagrams, 
so a lot is known about handling their quirks. 

2. (ANTI)FERROMAGNETS AND SPIN 
GLASSES 

For the Ising ferromagnet on three-regular (4> 3 ) 
graphs, solving the saddle point equations at large 

?! 
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(4) 



shows that the critical behaviour appears as an 
exchange of dominant saddle point solutions to 
the saddle point equations. The high and low 
temperature solutions respectively are 

= J+l 
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which give a low temperature magnetized phase. 
The critical exponents for the transition can also 
be calculated in this formalism and, as expected, 
are mean field. In general a mean field transition 
appears at 



exp(20 FM ) = z/(z - 2) 



(G) 



on cj) z graphs, which is the value predicted by the 
standard approaches. Simulations nicely confirm 
this mean field picture for the ferromagnet JTo| . 
Analysis of the Binder's cumulant for the magne- 
tization also shows that the critical temperatures 
are identical to the corresponding Bethe lattices 
(ie g = 3 for cf> 3 graphs). The specific heat is 
shown in Figure. 1 for various sizes of 4> 3 graphs. 
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Figure 1. The specific Heat for the Ising ferro- 
magnet on (j) 3 graphs of various size. 



There are various possibilities for addressing 
spin glass order in the Feynman diagram ap- 
proach. In H the entropy per spin was calculated 
for the Ising anti-ferromagnet on 4> 3 graphs and it 
was found to become negative for sufficiently neg- 
ative /3, which is often indicative of a spin glass 
transition. Simulations again confirm the picture. 
Taking a quenched distribution of couplings of the 
form 



P(J)=p5(J-l) + (l-p)6(J + l), 



(7) 



which gives the antiferromagnet for p = 0, pro- 
duces results for the spin glass order parameter, 
the overlap, that are very similar to the infi- 
nite range mean field (Sherrington-Kirkpatrick) 
model. Defining the overlap as 



(8) 



with two Ising replicas on each graph, <7j,Tj and 
histogramming 



P n (q)= (S(q-l/n 



i)> 



(9) 



where [ ] denotes the quenched disorder average, 
we get the distribution shown in Figure. 2 10 
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in the putative spin glass phase at low tempera- 
ture. The long tail stretching down to q = is 
characteristic of the mean-held spin glass picture 
of many inequivalent states. 
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Figure 2. P(q) at f3 = 1.2 for various p on </> 3 
graphs of size 500. 



For three or more replicas one does not see the 
continuous transition because a hrst order transi- 
tion occurs at higher temperature to a replica- 
symmetric state. The situation appears to be 
rather similar for Q > 2 state Potts models where 
the saddle point calculation finds a continuous 
transition at one of the spinodal points and misses 
a first order transition occurring at higher tem- 
perature. The 3-state Potts model, for instance, 
with action 



S = - 



+ (/> b + <j> c ) - c{4) a 4>b + (pa4> c + 4>b4> c ) 

(12) 



+ 4>i + <t>c), 



gives high and low temperature solutions 



>J a,b,c 



»a,b 



1 - 2c 



1 + VI - 4c - 4c 2 



1 + 2c - y/l-Ac- 4c 2 



(13) 



where c = 1/(^+1) 



It is possible to make some analytical inroads 
as well by looking at the solutions to the saddle 
point equations for k Ising replicas jll|] 



P(J) 



) k KdJd? 



(10) 



where we have denoted the 2 fc fields that now ap- 
pear as <j>. The Hessian for these equations is 
analytically calculable for any k 



n 



2 / P(J) tanh(/3J) m d/ - 1 



(A) 



(11) 



and its zeroes show that the k = transition tem- 
perature observed in simulations or calculated by 
using the analogy with the Bethe lattice is identi- 
cal to the k = 2 transition temperature. This 
also occurs in the finite replica version of the 
Sherrington-Kirkpatrick model p^| , so yet again 
the thin graph results are resolutely mean field. 




Figure 3. P(q) at = 2.2 for the 3-state Potts 
model for cj) 3 graphs of size 1000 (labelled) and 
500 (unlabelled). 
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Calculations and simulations for Potts glasses 
are just as easy as for the Ising spin glass. In 
Figure. 3 overleaf we have plotted the distribution 
of overlaps for the three state Potts model at low 
temperature. The overlap for a Q state Potts 
model is now defined as 



n 

71 ' 



(14) 



and the lack of an inversion symmetry in the 
spins gives a different pattern of replica symmetry 
breaking in mean-field theory to the Ising model. 
The numerical results are still consistent with a 
mean field picture. 

3. Conclusions 

In summary, spin models on thin graphs offer a 
promising arena for the application of ideas from 
matrix models, large-n calculations in field the- 
ory and bifurcation theory. In the spin glass case 
the tensor (or near-tensor) product structure of 
the inverse propagator allows some quite general 
expressions to be derived for the Hessian in the 
saddle point equations and offers a powerful line 
of attack on questions such as replica symmetry 
breaking. As a subject for numerical simulations 
they offer the great advantage of mean field re- 
sults with no infinite range interactions and no 
boundary problems. 

The bulk of the simulations were carried out 
on the Front Range Consortium's 208-node Intel 
Paragon located at NOAA/FSL in Boulder. CFB 
is supported by DOE under contract DE-FG02- 
91ER40672, by NSF Grand Challenge Applica- 
tions Group Grant ASC-9217394 and by NASA 
HPCC Group Grant NAG5-2218. CFB and 
DAJ were partially supported by NATO grant 
CRG910091. 
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